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A Method for Transforming a Digital Signal from the 
Time Domain to the Frequency Domain 

Background of the Invention 

Reversible transforms that map integers to integers are important tools for state-of-the-art 
lossless coding schemes. For integer inputs, such a transform generates integer outputs 
that approximate its floating-point prototype transform outputs. The original integer 
inputs can be completely recovered by passing the integer outputs from the forward 
transform through the inverse transform. 

A transform called integer modified discrete cosine transform (IntMDCT) is recently 
proposed and used in the ISO/IEC MPEG-4 audio compression. The IntMDCT can be 
derived from its prototype - the modified discrete cosines transform (MDCT). In the 
book [1], Malvar gave an efficient realization of MDCT by cascading a bank of Givens 
rotations with a DCT-IV block. It is well known that Givens rotation can be factorized 
into three lifting steps for mapping integers to integers [2]. Therefore, the realization of 
IntMDCT relies on an efficient implementation of integ^ DCT-IV. 

Integer transforms can be directly converted from their prototypes by replacing each 
Givens rotation with three lifting steps. Because in each lifting step there is one rounding 
operation, the total rounding number of an integer transform is three times the Givens 
rotation number of the prototype transform. For discrete trigonometric transforms 
(Discrete Fourier Transform (DFT), Discrete Cosine Transform (DCT), Discrete Hartley 
Transform (DHT), and Discrete W Transform (DWT)), the number of Givens rotations 
involved is usually at JVlogjiV^ level, where N is the block size. Accordingly, the total 

roimding number is also at iVlogjiV' level for the family of directly converted integer 

transforms. Because of the roundings, an integer transform only approximates its 
floating-point prototype. The approximation error increases with the number of 
roundings. 
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In this invention, a method for transforming a digital signal from the time domain to the 
frequency domain^ in particular reversible integer transform like DCT IV, DWT IV and 
other transforms. The total rounding number of the method according to the invention 
can be significantly reduced, for example in case of DCT IV to be as low as As a 
result, the approximation error of the method according to the invention is far less than 
that of the directly converted integer transforms. The computational complexity of the 
proposed method is also very low. 

Summary of the Invention 

The invention relates to a method for transforming a digital signal from the time domain 

« 

to the frequency domain. The method can be used for any types of digital signal, such as 
audio, image or video signals. The digital signal, which corresponds to a physical 
measured signal, may be generated by scanning at least a characteristic feature of a 
corresponding analog signal (for example, die luminance and chrominance values of a 
video signal, the amplitude of an analog soimd signal, or the analog sensing signal from a 
sensor). The digital signal comprises a plurality of data symbols. The data symbols of 
the digital signal are grouped into blocks, with each block having the same predefined 
number of data symbols based on the sampling rate of the corresponding analog signal. 

The digital signal is transformed to tiie frequency domain by a transforming element 
based on a transformation frmction comprising a transformation matrix. In this method 
according to the invention, each block of the digital signal is divided into two sub-blocks 
with equal number of data symbols. The transforming element transforms the digital 
signal block by block by processing, for each block, the two sub-blocks of data symbols 
as input signals simultaneously to generate two corresponding sub-blocks of transformed 
output signal, which can be combined to form a corresponding block of an output signal. 
It should be noted that the block size of the input signals and the output signals are the 
same. 

The transformmg element includes a plurality of lifting stages. A first lifting stage 
-rcceives-flie two sub-blocks as first and second input signals, respectively. The first input 
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signal is processed in a transformation path in the first lifting stage which includes a 
domain transformer, a rounding unit and a summation unit. The domain transformer can 
be any type of non-integer (real value) transformation function. The transformed first 
input signal is sununed with the second input signal and result in an ou^ut signal. The 
output signal and the first input signal of &e first lifting stage are used as the first and 
second input signals for a second (or subsequent) lifting stage, resulting in an output 
signal by the second lifting stage. Similarly, the output signal and the first input signal of 
the second lifting stage are used as the first and second input signals for another 
subsequent lifting stage. Since the output signal and flie first input signal fix)m each stage 
are received by the subsequent stage, these two signals shall be considered to be signals 
output from each of the lifting stages, although the second input signal of the subsequent 
stage is identical to the first input signal of ttie previous stage (i.e. not tratisformed). 

The transforming element can be illustrated based on the model of a lifting ladder. The 
lifting ladder model has two side pieces, each for receiving one of the two sub-blocks of 
data Sjmibols. Two or more cascading lifting stages are provided between the two side 
pieces. Each lifting stage receives a signal at one end (input end), and outputs a signal at 
the other end (output end) via a summation imit. A domain transformer is arranged at die 
input end of the lifting stage, and a rounding unit is arranged at the output end, between 
the domain transformer and the summation unit. The lifting stages are arranged between 
the side pieces in an alternating manner, such that the output (or input) ends of adjacent 
lifting stages are connected to the different side pieces. 

The number of lifting stages of the transforming element is defined by the number of 
lifting matrices which is determined by a process according to this invention as described 
in the following. 

This process first factorizes the transfomiation matrix of the transformation function into 
a string of first lifting matrices. Each lifting matrix comprises four sub-matrices, with 
two permutation matrices as two of the sub-matrices in one diagonal, and a zero and any 
general matrix as the other two of the sub-matrices in the other diagonal. A permutation 
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matrix is a matrix which changes the position of the elements in another matrix. An 
example of a lifting matrix is: 



'P2 



N 



0 



wherein is a general matrix, and Pl^ and P2j^ are permutation matrices. The 

permutation matrices may also be identical. Also, Identity Matrix is commonly used as a 
form of permutation matrix. An Identity Matrix is a matrix shown in the following: 



1 



1 



N 



1 



The general matrix Gj^ may take the form of any matrix configuration, for example any 
transformation matrix, with order 

Subsequently, the process identifies and combines any adjacent pairs of the fu^t lifting 
matrices in said string having their permutation matrices and the general matrix at the 
same sub-matrices position in each lifting matrix. Considering the following example: 



P2s 


0 • 






0 




A 






B 


P3„ 





P2„P4„ 



AP4„ +Pl„B PI 



, wherein A and B are 



0 

any general matrix. The two first lifting matrices (l^A side of the above equation) can be 
combined to form one second lifting matrix (right side of the equation), as the 
permutation matrices and the general matrix of both the first lifting matrices are in the 
same position. It should be noted that the product of two permutation matrices will also 
result in a permutation matrix. 
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As a result, a string of second lifting matrices is generated. It should be noted that the 
string of second lifting matrices may comprise a mixture of lifting matrices which are 
identical to the first lifting matrices and matrices which are a result of combination of 
first lifting matrices. The number of second lifting matrices then corresponds to the 
number of lifting stages in the lifting ladder of die transforming element. 

It should be noted that although the transformation element is described in the form of the 
lifting ladder model, it is only to illustrate the transformation paths of the transformation 
element However, the invention shall not be limited to said ladder model. 

Figure 1 shows a flow chart of an embodiment of the method according to the invention 
wherein the transformation matrix of the transformation function is converted by the 
above process into five second lifting matrices, resulting in the use of five lifting stages. 
From Fig.l, xj and X2 are first and second sub-blocks of a block of the digital signal, 
respectively, zj, zz and 25 are intermediate signals, and yi and y2 are output signals 
corresponding to the respective first and second sub-blocks. 

The method according to the invention can be used for transforming an input digital 
signal which represents integer values to an output signal which also represents integer 
values. In other words, a '^same word-length*' transformation is performed such that 
word-lengths of the iiq)ut signal and the output signal are the same (for example, a 8-bit 
input data is transformed to a 8-bit output data). Since both the input and output signals 
have tiie same word-length, the transformation method according to the invention is 
reversible. The output signal may be transformed back to the original input signal by 
performing the transformation mediod according to the invention. Such a reversibility 
property of the transformation according to the method of the invention can be used in 
lossless coding in which the output signal should be identical to the original input signal. 

Due to this reversibility property, the invention shall also include the transforming of a 
digital signal firom the frequency domain to the time domain using the transformation 
element as described earlier. 
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Such a same word-length (or integer) transfoimation of signals according to the invention 
can be used in many applications and systems such as MPEG audio, image and video 
compression, JPEG2000 or spectral analyzers (for analyzing Infrared, Ultra-violet or 
Nuclear Magnetic Radiation signals). It can also be easily implemented in hardware 
systems such as in a fixed-point Digital Signal Processor (DSP), without having to 
consider fectors such as overflow in the case of a real-value signal transformation. 

As can be seen from the method according to the invention, all die data symbols in each 
block of the digital signal are provided to the transforming element as a data vector. In 
each lifting stage, the data vector is transformed in the domain transformer, and the 
transformed data vector is rounded subsequently to an integer vector (i.e. after the 
transformation in each lifting stage). In other words, rounding is performed in the 
method of the invention once on the data vector as a whole. This is in contrast to any 
method according to the state of the art, wherein the rounding process is performed 
within the transformation process for the individual element or data symbol in each block 
of the digital signal. Thus, the number of rounding operations in the method according to 
the invention is greatly reduced. Due to tfie reduced number of rounding operations, the 
method according to the invention does not require large computation time and computer 
resource. 

Discrete Cosine Transforms. Discrete Sine Transforms Discrete Fourier Transforms or 
Discrete Wavelet Transforms are examples of transfoimation fimctions that may be used 
as the transformation fimction according to the invention. It should be noted tiiat the 
number of lifting stages of the transforming element may be different, depending on the 
result of tiie above described process carried out on the respective ti-ansformation 
functions. 

The factorizing of the transformation matrix into the string of die first lifting matiices 
may include decomposing the transformation matrix into a string of transformation sub- 
matrices which tiien in turn are decomposed into said string of first lifting matrices. 
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The corresponding transformation matrix of the transformation function is in general 
decomposed into three transformation sub-matrices. Such is the case, for sample when 
DCT-IV or DWT-IV is flie transformation function, wherein die detail implementation of 
both will be described in detail later. Each of the transformation sub-matrices are then 
further fectorized or decomposed into three lifting matrices, thereby forming nine lifting 
matrices as the string of first lifting matrices. As ah^ady described above, the nine lifting 
matrices are further combined, if possible, to form the string of second matrices. 
However, in some cases, less than nine lifting matrices forms the string of first lifting 
matrices when one or more of the transformation sub-matrices is a permutation matrix. 

For the case with DCT-IV as the transformation function, the three transformation sub- 
matrices corresponds to a pre-rotation matrix, a DCT-IV transformation-based matrix and 
a post-rotation matrix. In the special example, the DCT-IV transformation-based matrix 
and the post-rotation matrix can be combined to form a T matrix. Thus the factorization 
of the pre-rotation matiix and the T matrix only generate six lifting matrices as the stiing 
of first lifting matrices. As can be seen in one of the examples later, the last lifting matrix 
of the pre-rotation matrix is combined witii the first lifting matrix of the T matrix, 
resulting in a total of five lifting matrices in die string of second lifting matiices, the first 
two and the last two second lifting matrices are identical to those of the first lifting 
matrices. 

For the case widi DWT-IV as die tiansformation function, the three transformation sub- 
matrices corresponds to a rotation matrix, a DWT-IV transformation-based matrix and a 
permutation matrix. In this case, the permutation matrix is not fectorized into lifting 
matrices, and is implemented as a permutation block which rearranges the order of the 
elements of the input signal. Thus only the rotation matrix and the DWT-IV 
transformation-based matrix are factorized, generating six lifting matrices as tiie string of 
first lifting matrices. As can be seen again in one of the examples later, the last lifting 
matrix of die rotation matrix is combined with the first lifting matrix of the DWT-IV 
transformation-based matiix, resulting in a total of five lifting matiices in the stiing of 
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second lifting matrices, the first two and the last two second lifting matrices are identical 
to those of the first lifting matrices. 

It should be noted that each of the elements in the lifting matrices is a matrix itself. 

The invention not only relates to a method for transforming a digital signal fi-om the time 
domain to the firequency domain (and vice versa), but also includes a computer program, 
a computer readable medium and a device for implementing the said method 

Example of the method according to the invention based on an integer-type 
DCT IV as the transformation function. 

An example of the method according to the invention is illustrated using the DCT-IV 
transformation function. The DCT-IV of a //-point real input sequence jc(/i) is defined 
as follows [2]: 




(1) 



Let be the transformation matrix of DCT-IV, that is 



(m -t- l/2)(yi + V2>r 
N 




(2) 



The DCT-IV matrix can be factorized into the following form: 




R 



(3) 
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where i?^ is a x iV^ pre-rotation matrix defined by 



(4) 



i 

and /^/2 identity matrix of order 



Cl^,2 is the DCT-IV matrix of order N/2 and D^,2 is an order iV/2 diagonal matrix given 
by 




J?^ is a iSTx iST post-rotation matrix defined by 
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(6) 
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C"' D 



R- = 



1 



>/2 



C"' D 



-C"" D 



(7) 



Equation (3) is then repressed as 



(8) 



Both and 7* can be fiictorized into the product of three lifting matrices: 



fin 



Nil 



(9) 



where J5:,=-fez)^,,+V2/^,,):/r,, J^2-^ = VJc^", A/2 + . and 



po 





r/ 1 
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^Nfl ^ 2 


m 


^Nf2 








^NI2 

» m 




^3 ^Nn 



(10) 



where 



- tan 



- tan 



- tan 



JT 

SAT 
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Sin 



AN 



sin 



In 
4N 



sin 

4iV 



Thus, Equation (3) can be re-written as 



(11) 



Matrix S is defined as the product of and r, , that is 



•*JV/2 1 



(12) 



Obviously, S is also a lifting matrix. From (11) and (12), the final &ctorization fonnula 
for DCT-IV matrix is obtained as: 



(13) 



Equation (13) indicates that integer DCT-IV transform can be realized by five lifting 
stages. Because these five lifting stages have identical structures, only the transform 
equations for one stage is listed here. It is very easy and straightforward to write down 
the transform equations for all the other stages. The first stage 7^ is used as an example. 

Let * = [jc(«)Lo4... .AT-i J' " ('"^L-oj,. A-i ™P"* output Nxl integer 

vectors for the first transform stage: ' 



(14) 



P100254 



11 



vectors x and V are divided from the center into two halves, that is 







y. 








72. 



(15) 



Using (9) and (15), (14) is re-written as 



y, 








• 




y2. 




.^3 









(16) 



From (16), the time to frequency domain integer transform is: 



(17) 



where [*J denotes rounding operation. Figure 2 illustrates the method of transforming a 
digital signal from the time domain to the frequency domain using DCT-IV as the 
transformation function, xi and X2 are two sub-blocks of a block of the input digital 
signal, z/, 22 and zj are intermediate signals, and yi and y2 are corresponding sub-blocks 
of the output signal. 

Similarly, the frequency to time domain transform is given by: 



^2=y2-L^3-^ij 



(18) 



where [♦] denotes rounding operation. Figure 3 illustrates the method of transforming a 

digital signal from the frequency domain to the time domain using DCT-IV as the 
transfomiation function. 
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Example of the method according to the invention based on an Integer-type 
DOT IV as the transformation function. 

Another example of Ae method according to the invention is illustrated using the DWT- 
IV transformation function. The DWT-IV of a N-point real input sequence x(/i) is 
defined as follows: 



f2^V t \ • (ot + 1/2)(« + l/2)2jr \ 



1 



(19) 



Let Wl^ be the transformation matrix of DWT-IV, that is 



m,B-0,l,"-,/r-| 



(20) 



The DWT-IV matrix can be factorized into the following form: 



(21) 



is a N^cN rotation matrix 



R 



AT 



V2 



Nil 



Nil 



(22) 



In 12 is the identity matrix of order WLJ^^ is the counter identity matrix of order N/2 
which has the following form: 
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Pf, ISA NxN permutation matrix 



Nil 



N/2 



Tisa iV X JV^ matrix given by 



1 



^f2 



D 



IV 
N/2 



^N/2^N 



/2 



where CjJT/a is the DCT-IV matrix of order N/2 



^N/2 



Jn/2) 



cos 



( {m + l/2)(/i + l/2)jr V 



Dfii^ is an order N/2 diagonal matrix given by 



1 



-I 



D 



Nil 



1 



-1 



jtfoih R„ and Tcan be factorized into the product of three lifting matrices: 
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(28) 



where K, = (/2/^;2 - C/r/2i>Ar/2 )^«/2 ^2 = 



-C 



H/2 



V2 



r/ 




m 


^N/2 ^ 2 


r/ 1 




m 






^3 



(29) 



where 



- 0.414 



-0.414 



- 0.414 



N/2 



0.707 



0.707 



0.707 



Thus, Equation (21) can be re-written as 



(30) 



Matrix 5* as the product of and T^ , that is defined as: 



s = i?,?; 



r r 



(31) 
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Obviously, S is also a lifting matrix. From (30) and (31), the final fsLctorization formula 
for DWT-IV matrix is obtained as: 



(32) 



Equation (32) indicates that integer DWT-rV transform can be realized by five lifting 
stages plus a data shuffling stage defined by permutation matrix . The data shufiCling 

stage rearranges the components order in each input data blocic According to , the 

input data vector is rearranged in the following way: the first half of vector remains 
unchanged; the second half of vector is flipped — tiiat is 









• 




• 
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m 




• 








• 
• 
• 










• 
• 

• 









(33) 



Because the five lifting stages have identical structures, only the transform equations for 
one stage is listed here. It is very easy and straightforward to write down the transform 
equations for all the other stages. The first stage is used as an example. 

Let >f =[jc(«)L.oj,...jv-, and y = bC'n)]^.^^^..^., be the input and output Nxl integer 
vectors for the first transform stage: 



(34) 



Vectors X and y are divided from the center into two halves, that is 
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y, 








^2. 



(35) 



Using (28) and (34), (33) is re-vnritten as 



yi 




r/ 1 


• 















(36) 



From (36), the time to frequency domain integer transform is: 



(37) 



Figure 4 illustrates the method of transforming a digital signal from the time domain to 
the frequency domain using DCT-IV as the transformation function, jc/ and X2 are two 
sub-blocks of a block of the input digital signal, Zi, Z2 and zj are intermediate signals, and 
yi and>^2 are corresponding sub-blocks of the ou^ut signal. 

Similarly, the frequency to time domain integer transform is given by 



(38) 



Figure 5 illustrates the method of transforming a digital signal from the frequency 
domain to the time domain using DCT-IV as the transformation fimction 



Conclusion of the two examples 

From (17) and (18), there are N/2 roundings in each lifting stage. Therefore, from (13), 
the total rounding number for the proposed integer DCT-IV algorithm is five times N/2, 
that is 2.5N^ which is significantly lower than NJog^N according to the state of the art In 
(13), the majority of computation power is used in the four N/2 point DCT-IV 
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subroutines C^,2 » when N is a, large value, e.g. N = J 024, Because the floating-point 
DCT-IV Cj^ can be calculated using two half-length DCT-IV CjJ^^j plus pre- and post- 
rotations according to (3), the arithmetic complexity of the proposed integer DCT-IV is 
roughly estimated to be twice that of the floating-point DCT-IV. 

A similar conclusion can be drawn for the integer DWT-IV transformation function, 
wherein Ave lifting stages are used. 

In this invention, a method for realizing reversible transformation function, for example 
for integer type-IV DCT and DWT transformation function, is proposed. This method 
requires significantly reduced number of roundings for every block of input samples. 
As a result, the approximation error is greatly reduced. The method according to the 
invention is low in computational complexity and modular in structure. 
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What is claimed is: 

1 . A method for transforming a digital signal from ttie time domain into tide frequency 
domain using a transfomiation frmction comprising a transformation matrix, the digital 
signal is divided into a plmality of blocks, each block comprising a predefined number of 
data symbols, the method comprising 

- dividing each block into two sub-blocks with equal number of data symbols, 

- transforming the two sub-blocks by a transforming element, the transforming element 
comprises a plurality of lifting stages, wherein the two sub-blocks are received by a first 
lifting stage of the plurality of lifting stages as first and second input signals, and wherein 
each lifting stage comprises a transformation path, wherein the first input signal received 
by the lifting stage is processed by a domain transformer and a rounding unit, and is 
subsequently summed with the second input signal in a summation unit to result in an 
ou^ut signal; and wherein each sul^sequent lifting stage receives the output signal of the 
previous lifting stage as the first input signal^ and the first input signal of the previous 
lifting stage as the second input signal, to generate the corresponding output signal of the 
subsequent stage, and wherein the number of lifting stages is defined by the number of 
lifting matrices detemiined by a process comprising 

ftictorizing the transformation matrix into a string of first lifting matrices, each of 
which comprises two permutation sub-matrices in one diagonal, and a zero and a general 
sub-matrix in the other diagonal, 

combining all such pairs of first lifting matrices being adjacent in the said string 
the respective permutation sub-matrices of which are identical and the general sub- 
matrices of which are at their corresponding positions, thereby generating second lifting 
matrices, wherein the number of stages corresponds to the number of the second lifting 
matrices. 
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2. The method of claim 1, wherein the traiisfomiation function is a DCT I transfomiation 
function, DCT II transfomiation function, DCT IV transformation function, FFT 
transformation function or D WT-IV transformation function. 

3. The method of claim 1, wherein the step of &ctorizing the transformation matrix 
comprises decomposing the transformation matrix into a string of transformation sub- 
matrices and decomposing said transformation sub-matrices, if possible, into the said first 
lifting matrices. 

4. The method of claim 2, wherein the step of factorizing the transformation matrix 
comprises decomposing tiie transformation matrix into a string of transformation sub- 

« 

matrices and decomposing said transformation sub-matrices, if possible, into the said first 
liftmg matrices. 

5. The mediod of claim 3 or 4, wherein each of said decomposed matrices are furtiier 
decomposed into three lifting matrices, resulting in nine first lifting matrices. 

6. The method of claim 4, wherein the DCT IV matrix of the DCT IV transformation is 
decomposed into a pre-rotation matrix^ a DCT IV transformation-based matrix and a post 
rotation matrix as transformation sub-matrices. 

7. The method of claim 6, wherein the DCT IV transformation-based matrix and the post 
rotation matrix are combined into a T matrix, and the pre-rotation matrix and the T matrix 
are each decomposed into three first lifting matrices, resulting in six first lifting matrices. 

8. The method of claim 7, wherein the last lifting matrix of the pre-rotation matrix is 
combined with the first lifting matrix of the T matrix, resulting in a total of five second 
lifting matrices. 

9. A method according to any one of the preceding claims, wherein each element of the 
first and second lifting matrices is a matdx. 
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10. A device for transforming a digital signal from the time domain into the 
frequency domain using a transforaiation frmction comprising a transformation matrix, 
the digital signal is divided into a plurality of blocks, each block comprising a predefined 
number of data symbols, the device comprising 

- a dividing unit for dividing each block into two sub-blocks witti equal number of data 
symbols, 

■ 

- a transformation unit for transfomiing the two sub-blocks by a transforming element, 
the transforming element comprises a plurality of lifting stages, wherein the two sub- 
blocks are received by a first lifting stage of the plurality of lifting stages as first and 
second input signals, and wherein each lifting stage comprises a transformation path, 
wherein the first input signal received by the lifting stage is processed by a domain 
transformer and a rounding unit, and is subsequently summed with the second input 
signal in a summation unit to result in an output signal; and wherein each subsequent 
lifting stage receives the output signal of the previous lifting stage as the first input 
signal, and the first input signal of the previous lifting stage as the second input signal, to 
generate the corresponding output signal of the subsequent stage, and wherein the 
number of lifting stages is defined by the number of lifting matrices determined by a 
process comprising 

factorizing the transformation matrix into a sting of first lifting matrices, each of 
which comprises two permutation sub-matrices in one diagonal, and a zero and a general 
sub-matrix in the other diagonal, 

combining all such pairs of first lifting matrices being adjacent in the said string 
the respective permutation sub-matrices of which are identical and the general sub- 
matrices of which are at their corresponding positions, thereby generating second lifting 
matrices, wherein the number of stages corresponds to the number of the second lifting 
^matrices. 
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1 1 . A computer readable medium, having a program recorded thereon, wherein the 
program is to make the computer execute a procedure for transforming a digital signal 
from the time domain into the frequency domain using a transformation fimction 
comprising a transformation matrix, the digital signal is divided into a plurality of blocks, 
each block comprising a predefined number of data symbols, the method comprising 

- dividing each block into two sub-blocks witti equal number of data symbols, 

- transforming the two sub-blocks by a transforming element, flie transforming element 
comprises a pluraUty of lifting stages, wherein the two sub-blocks are received by a first 
lifting stage of the pluraUty of lifting stages as first and second input signals, and wherein 
each lifting stage comprises a transformation path, wherein the first input signal received 
by the lifting stage is processed by a domain transformer and a rounding unit, and is 
subsequently summed with Has second ii^mt signal in a summation unit to result in an 
output signal; and wherem each subsequent lifting stage receives the output signal of the 
previous Ufting stage as the first input signal, and the first input signal of the previous 
lifting stage as the second input signal, to generate the corresponding output sig^ of the 
subsequent stage, and wherein the number of lifting stages is defined by the number of 
lifting matrices determined by a process comprising 

ftictorizing the transformation matrix into a sting of first lifting matrices, each of 
which comprises two permutation sub-matrices in one diagonal, and a zero and a general 

sub-matrix in the odier diagonal, 

combining all such pairs of first lifting matrices being adjacent in the said string 
tiie respective permutation sub-matrices of which are identical and the general sub- 
matrices of which are at their conresponding positions, thereby generating second lifting 
matrices, wherein the number of stages corresponds to tiie number of the second Ufting 
matrices. ' 

12. A computer program element which is to make the computer execute a procedure for 
transforming a digital signal &om tihe timf! domain into the firequency domain using a 
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transformation function comprising a transformation matrix, the digital signal is divided 
into a plurality of blocks, each block comprising a predefined niunber of data symbols, 
the method comprising 

- dividing each block into two sub-blocks with equal number of data symbols, 

- transforming the two sub-blocks by a transforming element, the transforming element 
comprises a plurality of lifting stages, wherein the two sub-blocks are received by a first 
lifting stage of the plurality of lifting stages as first and second input signals, and wherein 
each lifting stage comprises a transformation path, wherein the first input signal received 
by the lifting stage is processed by a domain transformer and a roimding unit, and is 
subsequently summed with the second input signal in a siunmation unit to result in an 
ou4)ut signal; and wherein each subsequent lifting stage receives the output signal of the 
previotis lifting stage as the first input signal, and the first ii^>ut signal of the previous 
lifting stage as the second input signal, to generate the corresponding output signal of the 
subsequent stage, and wherein the niunb^ of lifting stages is defined by the niunber of 
lifting matrices determined by a process comprising 

factorizing the transformation matrix into a sting of first lifting matrices, each of 
which comprises two permutation sub-matrices in one diagonal, and a zero and a general 
sub-matrix in the other diagonal, 

combining all such pairs of first lifting matrices being adjacent in the said string 
the respective permutation sub-matrices of which are identical and the general sub- 
matrices of which are at their corresponding positions, thereby generating second lifting 
matrices, wherein the number of stages corresponds to the number of the second lifting 
matrices. 

r 
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1st lifting stage 



2nd lifting stage 
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Figure 1. 
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Figure 2. Forward IntDCT-IV flowchart 
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Figure 3. Inverse IntDCT-IV flowchart 
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Figure 4. Forward IntDWT-IV flowchart 
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Figure 5. Inverse InlDWT-IV flowchart 
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